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The magnetic susceptibility and susceptibility anisotropy of the quasi-binary alloy system 
Cea;Lai_a;Cu2.05Si2 have been studied for low concentration of Ce ions. The single-ion description 
is found to be valid for x < 0.09. The experimental results are discussed in terms of the degenerate 
Coqblin-Schrieffer model with a crystalline electric field splitting A ~ 330 K. The properties of 
the model, obtained by combining the lowest-order scaling and the perturbation theory, provide a 
satisfactory description of the experimental data down to 30 K. The experimental results between 
20 K and 2 K are explained by the exact solution of the Kondo model for an effective doublet. 



PACS 72.15.Qm, 75.20.Hr, 75.30.Mb 
A. Introduction 

The intermetallic compound Ce :r Lai_ :r Cu2.o5Si2 bas 
been studied for more than two decades, but its prop- 
erties are still not completely understood. For a 
large concentration of Ce ions ţhc compound shows 
heavy fermion supiepconductivitya and small-moment 
antifcrromagnetismflfl At intermediate concentrations, 
non-Fermi-liquid features are foundja while for small x 
one observes anomalies due to the crysiaJline electric field 
splitting (CF) and the Kondo effect.ErO The change of 
the properties of the compound induced by Ce doping 
is an important issue that remains to be clarified, both 
experimcntally and theoretically. 

Here we present magnetic susceptibility data of 
Ce K Lai_ x Cu2.o5Si2 intermetallics for small concentra- 
tions of Ce ions, x < 0.09, and show that the mag- 
netic mome&t of Ce depends on temperature and is very 
anisotropic.B The samples we consider are sufhciently di- 
lute for the interaction between the Ce ions to be ne- 
glcctcd. We find the susceptibility of a single Ce ion 
embedded in a tetragonal host by studying the effects 
of Ce concentration. The magnetic anomalies are ac- 
companied by transport anomalies, as indicated by large 
peaks in the thermoelectric power and a logarithmic in- 
crease of the electric resistance.Q This behavior is typical 
of exchange scattering on CF-split ions and we base the 
theoretical analysis on the degenerate Coqblin-Schrieffer 
modelEll in which the two lowest CF states are coupled 
by exchange interaction to conduction states. 

The CF parameters of dilute Ce alloys are obtained 
by analyzing the high-temperature magnetization data 
with the usual CF theory which neglects the coupling 
between the / electrons and the conduction states, and 
considers just an isolated Ce ion in a tetragonal envi- 



ronment (for details see the Appendix). The magnetic 
anisotropy is explained by taking A ~ 330 K for the CF 
splitting between the doublet ground state and a pseudo 
quartet, formed by the two closely-spaced excited dou- 
blets, and taking r\ = 0.816 for the relative weight of the 
CF orbitals. These values are consistent with the high- 
temperatup. neutron scattering data on CeCu2Si2 single 
crystals.EjO However, the CF theory does not describe 
correctly the response of the high-temperature local mo- 
ment, and it also fails to explain the rapid reduction of 
the local moment at intermediate temperatures and the 
loss of anisotropy at low temperatures. These charac- 
teristic features of the magnetic response of Ce ions in 
metallic hosts can only be obtained by coupling the / 
states to the conduction bând. 

In this paper we show that the magnetic properties of 
dilute Ce x Lai_^Cu2.o5Si2 compounds can be explained 
by describing the interaction between the Ce ions and 
the metallic electrons by the Coqblin-Schrieffer model 
with the CF splitting. The high-temperature suscepti- 
bility obtained by the lowest-order perturbation theory 
agrees very well with the experimental data, provided we 
renormalize-the exchange interaction by the 'poor man's 
scaling'. EjTLZI For large CF splitting, the scaling theory 
generates two relevant low-energy scales, Tr- -C T§, 
instead of one without splitting. The behavior of the 
anisotropic susceptibility above 30 K is accounted for 
by choosing T K = 8.5 K and Tg = 100 K. The two 
low-energy scales, which differ by an order of magnitude, 
explain the reduction of the local moment from a high- 
temperature sextet to a low-temperature doublet 
lead to a qualitative explanation of the transport data 
The low temperature properties of the Coqblin-Schrieffer 
model can not be calculated by scaling. However, for 
A = 330 K and Tg = 100 K, the occupation of the ex- 



cited CF states below 20 K is negligibly small. Thus, 
we approximate the / state by an effective doublet, and 
describe the low-temperature properties by an effective 
spin-1/2 Kondo model with Kondo scale Tk- The ep. 
act renormalization group results for the susceptibilitytij 
agree very well with the experimental data. 

The paper is organized as follows. First we describe the 
sample preparation and provide the details of the mea- 
surements. Then we analyze the susceptibility data ob- 
tained by the Faraday magnetometer on powdered sam- 
ples and the anisotropy data obtained by the torque mag- 
netometer on polycrystalline samples, and show that the 
observed anomalies can not be explained by the CF ef- 
fccts due to an isolated / state. The coupling between 
the / state and the conduction bând, as described by the 
Coqblin-Schrieffer model with a CF splitting, is intro- 
duced next. Then, we discuss the susceptibility and the 
magnetic anisotropy of an excited quartet separated from 
the ground state doublet by the energy A. Finally, com- 
pare the theoretical results the experimental data. The 
CF calculations for Ce ions in a tetragonal environment 
are presented in the Appendix. 



B. Experimental details and data analysis 



plotted as a function of temperature for two pairs of typ- 
ical samples. The susceptibility difference for the two 
samples with the same nominal concentration x — 0.07 
grows from about 2 x IO -6 (emu/ mol) at 300 K to about 
200 x 10~ 6 (emu/mol) at 2 K, while x{ x ,T) increases in 
the same temperature interval only about 10 times. Sim- 
ilar behavior is seen in othcr samples, and we take that 
as an evidence that the susceptibility difference shown 
in Fig. 1 is not due to an inhomogeneous distribution 
of Ce ions. The slight variation of the funcţional form of 
x(xjT) is most pronounced at low temperatures (see also 
Fig.ra) , and that can be explained in terms of a small vari- 
ation in the Kondo temperature Tk- The inset in Fig. 
1 shows the data for the two samples with nominally 7 
at% of Ce, plotted versus T/Tk- The low-temperature 
data follow the same curve, provided we use Tk = 7.6 
K and Tk = 9.3 K for the data pepresented by the open 
and filled triangles, respectivelyEil However, at interme- 
diate temperatures, the data represented by the open and 
filled symbols can not be represented by a single function 
of T /Tk- The Kondo temperature of various samples 
might differ because of small deviations of the actual Cu 
concentration from the nominal one and the associated 
chemical pressure efiects. 



1. The samples 

The polycrystalline Ce a; Lai_ 2 ;Cu2.o5Si2 samples with 
low Ce-content (x < 0.09) were prepared using a two- 
step procedure in order to enhance composition control 
and improve homogeneity. In the first step, appropriatc 
amounts of pure elements were arc-melted to get maş- 
ter ingots with compositions x = and x = 0.1. In 
the second step, part of these maşter ingots were mcltcd 
together in an appropriate ratio to obtain samples with 
0.01 < x < 0.09. Samples with larger Ce-content were 
prepared in a single step directly from pure elements. 
AII samples were annealed for 40 hrs at 700 °C, then 
for 80 hrs at 950 °C, followed by a slow cooling to 700 
°C within 72 hrs. The X-ray powder diffraction patterns 
showed only reflections of the ThCi'2Si2 structure. Both 
lattice parameters a and c, as well as the unit cell vol- 
ume, decrease linearly with increasing Ce content, the 
decrease being quite pronounced for a but rather weak 
for c.Q This preparation procedure combines all the ex- 
perience gained in our laboratory within the last twenty 
years upon the investigation of CeCu2Si2 and related al- 
loys. More detailed information about the prepacation 
process can be obtained directly from the authors.£j 

The temperature dependence of the susceptibility 
x(x, T) is found by measuring the fixed powders weigh- 
ing about 20 mg. The data analysis shows that x( x i T) is 
dominated by the contribution due to Ce ions, but that. 
the samples with the same nominal Ce ConcentrationE2l 
show slightly different susceptibilities at low tempera- 
tures. This is demonstrated in Fig. 1, where x(x,T) is 
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FIG. 1. x( x jT) plotted versus temperature for two pairs 
of Ce 2! Lai- :c Cu2.o5Si2 samples. The triangles correspond to 
x = 0.07 and the squares toi = 0.03, respectively. The inset 
shows x( x j T)-Tk for the two x — 0.07 samples, plotted versus 
T/Tk , with two different values of Tk ■ 

The samples used for the torque measurements are 
small polycrystallites of irregular shape. The uniaxial 
symmetry of the unit cell allows one to obtain the intrin- 
sic anisotropy even in the absence of perfect alignmcnt 
of the crystallites (for details see below). 



2. The susceptibility 

The absolute value of the average magnetic suscepti- 
bility x{ x i T) is measured for Ce a; Lai_ a ;Cu2.o5Si2 samples 
on a high-precision Faraday balance for tempcratures be- 
tween 2 K and 350 K, for the magnetic field up to 9 
kOe, and for Ce conccntration ranging from 1 to 9 at%. 
The magnetization of all samples, measured at 300 K, 
77 K, and 2 K, is found to be a linear function of the 
applicd field. Susceptibility is defined as x( x , T) = M/H 
[emu/mol]. The experimental rcsults arc plotted in Fig. 
2 as a function of tempcrature and in Fig. 3 as a function 
of concentration. Fig. 3 shows that the susceptibility in- 
crease linearly with Ce concentration, and we take that 
as an evidence that the Ce- Ce intcraction is small. 
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FIG. 2. Susceptibility of Cea;Lai_a;Cu2.05Si2 plotted versus 
temperature on a logarithmic scale for x < 0.09. 

The data analysis shows that the high-temperaturc 
susceptibility of the samples with the lowest concentra- 
tion of Ce is diamagnetic, i.e., the response of very di- 
lute alloys is dominated by a background. Thus, the 
cffccts of the alloying on the matrix in which the Ce ions 
fluctuate should not be a priory neglected. We also no- 
tice that the low-temperature susceptibility has a Curic- 
like upturn due to some unspecified magnetic impuritics 
that are immancnt to the rare earths, and for very di- 
lute alloys this upturn becomes relatively large. Because 
the concentration of these unspecified magnetic impu- 
ritics diffcrs in cach samplc, and becausc the alloying 
changes the background, we do not define the single-ion 
Ce-response as a difference between the susceptibility of 
Ce 2 ,Lai_ a; Cu2.o5Si2 and LaCu2.osSi2. To minimize the 
systematic error due to sample preparation we define the 
susceptibility of a Ce-free compound as a statistical av- 
erage obtaincd from Fig. 3, and find the single-ion con- 
tribution by the following procedurc. 
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FIG. 3. Susceptibility of Ce x Lai_ a ,Cu2.o5Si2 plotted versus 
concentration for various temperatures. 

We assume that the total susceptibility of a given 
samplc, which is shown in Figs. 2 and 3, comprises 
two contributions: a single-ion Ce susceptibility, Xion, 
and a concentration-dependent background, xb(x,T) ~ 
%(0, T) + x ■ [8xb{x) I 'dx] x=0 . We consider only the low- 
est order concentration effects and express the measured 
susceptibility of a given sample as, 

X (x, T) - X (0, T) + x { Xt on(T) + [d X B(x)/dx] x=0 } . (1) 

To find xb(x,T) we notice that the susceptibility data 
shown in Fig. 3 are linear in concentration for all tempcr- 
atures, and use the x = intercept of the straight lines to 
define the statistically averaged background susceptibil- 
ity of the Ce-free matrix, ^(Ojî 1 ). The values obtaincd 
in such a way are shown in Fig. 4, together with the 
lcast-squarcs fit based on the expression 



X(0,T)= X o + 



a 



trup 



T-e t 



(2) 



where xo describes the constant diamagnetic contribu- 
tion of the stoichiometric compound LaCu2.05Si2 and 
the Curie- Weiss term describes magnetic contamination. 
The experiment gives, xo = —3.1 • 10~ 5 (emu/mol), 
Cimp — 5.9 • 10 -4 K (emu/mol), and Qi mp = —2.8 
K. To estimate the level of contamination we assume 
Cimp = {N lmp /3k B )g 2 fj, 2 B j(j + 1), where \x B is Bohr 
magneton, g is the gyromagnctic factor and kg is Boltz- 
mann constant. Taking for j the smallest possible value, 
j = 1/2, we find that the upper limit for the number of 
magnetic impurity atoms is Ni mp = 0.0016 Na, where 
Na is Avogadro number. Thus, the average concentra- 
tion of the unspecified magnetic impurities is about 0.1 
at%, which is at least one order of magnitude less than 
the lowest Ce concentration we measurc. 




Temperature (K) 

FIG. 4. The susceptibility of the LaCu2.o5Si2 matrix 
x(0, T) plotted versus temperature (circles). The full line 
is the least-squares fit by Eq. (0). The low-temperature 
Curie- Weiss upturn is due to magnetic contamination. Note 
the difference in the vertical scales used here and in Fig. 2. 

Using Eqs. (0) and (|) we find that [x(x, T)~x(0, T)]/x 
defines the same curve for all the samples. This uniquc 
curve defines the impurity susceptibility Xion(T) up to a 
constant shift [dxB{x)/dx] x _ . However, the two terms 
in the braces in Eq. (|lj) can not be obtained by inde- 
pendent measurements, and we replace Xion(T) in Eq. 
([j]) with a theoretical expression (for details see below) 
and find that the experimental data can be fitted by tak- 
ing [dxB/dx] x=0 = 10~ 4 (emu/mol Ce). The same value 
of this constant shift is obtained by using for the high- 
temperature single-ion susceptibility in Eq. (pi) the ex- 
pression Xion(T) ~ l/T, and making the least-squares 
fit through the data above 250 K. The comparison of 
XioniT) and dxB{x)/dx shows that the susceptibility 
contribution due to Ce ions is at least one order of magni- 
tude larger than the background contribution, so that a 
possible small error in [9xb(x) / dx] x _ does not influence 
our conclusions in an essential way. 

The average susceptibility of a single Ce ion in dilute 
Ce K Lai_ x Cu2.o5Si2 samples is defined by the expression 



Xion\l ) 



X (x,T) - X (0,T) 



IQ~ 4 (emu/ mol Ce), (3) 



which is shown in FigJq as a function of temperature for 
various values of x. The response of Ce ions is nearly the 
same for all the compounds up to x ~ 0.09, and it begins 
to deviate from the cc-independent form for x > 0.1. u The 
scatter of the data at low temperatures is about the same 
as the scatter shown in Fig. 1, and can be eliminated by 
plotting Xion on a universal temperature scale T/Tk(x), 
with Tk(x) — 8.5 ii K. However, in our dilute sam- 
ples, the lack of any systematic behavior of Tk{x) as a 
function of Ce-doping does not allow us to explain the 



observed fluctuations in terms of the chemical pressure ef- 
fects induced by Ce doping. As discussed in the previous 
section, the random variation of Tk{x) can be associated 
with the local fluctuations in Cu stoichiometry. 

The anomalous behavior of the molar susceptibility of 
Ce in Ce x Lai_ a; Cu2.o5Si2 becomes most transparent if 
one compares Xion(T) with the susceptibility of a crys- 
talographically equivalent but unhybridized 4/ 1 state, 
Xcf(T), which is calculated in the Appendix. In FigJŞjwe 
plot the susceptibilities versus logT, and in Fig.||we use 
the Curie- Weiss plot and show the inverse of the suscep- 
tibilities as a function of T. The symbols represent the 
experimental data, whilc the dashed line, the full line and 
the dashed-dotted line represent the CF, the scaling, and 
the Wilson's result, respectively. 
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FIG. 5. The av- 

erage single-ion susceptibility of Ce :r Lai_ :r Cu2.05Si2 defined 
as X ion(T) = [x(x,T) - x(0,T)]/x - 1 ■ lQ- A (emu/mol Ce), 
plotted versus temperature for x < 0.09 samples. The full 
line shows the scaling result, the dashed line shows the CF 
result and the dashed-dotted line shows the exact solution for 
the spin- 1/2 Kondo model in the local moment regime. 

The high-temperature susceptibility follows the Curie- 
Weiss law, l/xi<m(T) ~ (T - 9)/C, where C = 
NA(i% s /3kB — 0.125[/% s /[/,% ~ 0.8 K emu/mol Ce, which 
is close to the CF result. The antiferromagnetic Weiss 
temperature is ~ —100 K. The low-temperature sus- 
ceptibility data can also be represented by a Curie- Weiss 
law with C — 0.3 K emu/mol Ce and ~ — 15 K. This 
value of C leads to /x c ff ~ 1.6/is which is close to the av- 
erage value obtained for the lowest CF doublet (see the 
Appendix) . 

The analysis of the single-impurity data in Figs. || and 
o shows clearly that the CF theory, which neglects the 
coupling of / electrons to conduction states, fails to ex- 
plain the response of the 4/ 1 ions to an externai magnetic 
field. The correct description is obtained by considering 



not only the thermodynamic fluctuations but taking also 
into account the quantum fluctuations, due to Kondo 
effect. Before evaluating these effects in more detail, we 
discuss the torque measurements which provide the mag- 
netic anisotropy of Ce ions and allow us to obtain the full 
susceptibility tensor. 
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FIG. 6. The inverse of the average single-ion susceptibility 
shown as a function of temperature for Ce ;c Lai_; C Cu2.o5Si2 
samples with x < 0.09. Dashed line is the CF result, and the 
full line is the scaling result. 



3. The susceptibility anisotropy 

The in- 

trinsic anisotropy of polycrystalline Ce a: Lai_ 2: Cu2.o5Si2 
samples with clear preferenţial orientation was measured 
on a high-prccision torque magnetomctcr. In the exper- 
iment the sample is attached to a thin vertical quartz 
fiber and a homogeneous magnetic ficld H is applied in 
the horizontal xy plane, as shown in Fig. 7. Thus, only 
the z component of the torque T is measured. 



c-* 




FIG. 7. The torque on a uniaxial crystallite, oriented along 
c direction, induced by the externai magnetic field H in the 
xy plane. 

The experimental values for the susceptibility 
anisotropy of a given sample Ax exp are obtained from the 
ratio T z /H 2 , which is measured by the following proce- 
dure. For an arbitrary sample orientation, the magnetic 
field is rotated in the xy plane by an anglc a varying 
from 0° to 180°. Measuring the angular dependence of 
the torque curve at 300 K, 77 K and 2 K, in a mag- 
netic field of 8 kOe, we find that T z (a) is a sinusoidal 
function of a period 7r, with zeros which are indepen- 
dent of temperature. This shows that the magnetiza- 
tion of the sample is proporţional to the applied field. 
The amplitude T z of the sine torque curve is then mea- 
sured as a function of temperature. Next, the sample 
is rotated by 90° around the y axis, and the amplitude 
T z v is measured as a function of temperature. The ex- 
perimental anisotropy of a given sample is defined as 
^xt V = 2IT z I H 2 n [emu /mol], where n is the number 
of moles and the index i — h, v denotes the two orien- 
tations of the sample. The data obtained in this way 
are shown in Fig. 8 where Ax^ p (x,T) is plotted versus 
temperature for x = 0.01, 0.02, 0.03, and 0.06. AU of the 
curves in Fig. 8 show similar temperature dependence, 
with the extremum point at 48 K. The torque signal is 
quite strong due to a large anisotropy of the Ce moment 
and because the measurements are performed on samples 
with a rather high degree of preferenţial alignment of the 
crystallites within the polycrystal. 
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FIG. 8. The susceptibility 

anisotropy of Ce ;c Lai_ :E Cu2.05Si2 is shown as a function of 
temperature for x < 0.06 samples. The data corresponding 
to two difîerent orientations of the sample with respect to the 
magnetic field are denoted by h and v, respectively. 



To analyze the torque data shown in Fig. 8 we write 



the z component of T as 



T z = V(M x H) 



(4) 



where V is the volume of the sample and M is the mag- 
netization induced by H. For a single tetragonal crystal- 
lite, the induced magnetization has an isotropic and an 
anisotropic component, i.e., it can be written as 



M = M J 



•M^ 



(5) 



where the isotropic component M 7 is directed along H, 
while the anisotropic one M is directed along the c axis 
of the crystal. Obviously, M 7 = x ab H and M A = (x c - 
X a6 )H c , where x c 1S the susceptibility along the c axis, 
X ab is the susceptibility in the ab plane, and H c is the 
component of the magnetic field along the c direction. 
Using Eq.([|), we find that the torque on a crystallite 
with the c axis directed along (9, (p), and induced by the 
magnetic field in the xy plane, H = H(cosa, şina, 0), is 
given by the expression (see also Fig. 7) 

T z = ^nH 2 A X sin 2 (9) sin [2(a-ip)], (6) 

where A% is now measured in emu/mol, i.e., stands for 
(V/n)Ax- Thus, the temperature dependence of the 
torque is given, up to a numerical constant, by the intrin- 
sic susceptibility anisotropy Ax = X° ~ X ah , regardless of 
the orientation of the crystallite. 

In a polycrystalline sample, we are dealing with somc 
volume distribution Vi = Vi(8i, (fi) of the crystallites with 
respect to their orientations. The total torque on a sam- 
ple with a given concentration of Ce ions is the sum of 
all the separate contributions and can be written as 



K = ~nH 2 lx A X sin[2(a-<p x )}, 



(7) 



where we introduced the alignmcnt factor j x (0 < j x < 
1) and the phase shift (p x . In our experiment the magnet 
is rotated to a position a — Cp x — ±tt/4: which maximizes 
the torque. Thus, up to a geometrical factor j x , the 
temperature dependence of the torque follows from the 
intrinsic anisotropy of the sample. The alignment factor 
7a; depends on the distribution of the polycrystallites in 
the sample, and is given by j x — y/A 2 + B 2 , where A x = 
yJ2%Vi sin2 l cos2(p i and B x = ± Y^i v i sin2 ^sin2^j. 
To find the single-ion anisotropy we have to estimate ^ x 
and the intrinsic anisotropy of the metallic matrix. 

The data analysis shows that the experimental 
anisotropy Ax exp (x,T) = 7a;Ax follows a similar tem- 
perature dependence in all the samples. This is shown 
in Fig. 9, where Axi Xp (x,T) is plotted as a function of 
Axy Xp {xi,T), with temperature as an implicit variable. 
We recall that index i — v,h denotes the orientation of 
the sample with respect to the applied field. The orien- 
tation of the reference sample is v and the referent Ce 
concentration is x\ = 0.01. The correlation between the 
torque data is found to be strictly linear, 



AxŢ v {x,T)=kiA X l*r{x x ,T) + li 



(8) 



where k x and l x denote the slope and the intercept of the 
straight lines respectively. The numerical values of k x 
and l x for each sample are given in Table I. As shown 
below, the linear relation (@) proves essential for ob- 
taining the single-ion contribution from the experimental 
anisotropy data. 
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FIG. 9. Axl xp (x, T) of Ce a; Lai_; C Cu2.o5Si2 samples is plot- 
ted versus Axv Xp ( x i,T). The data are given for two different 
orientations (i = h, v) with respect to H. The concentration 
of the Ce ions in the referent sample is xi — 0.01. The full 
lines are the linear extrapolations which define the coefficients 

In analogy with the analysis of the susceptibility data 
we assume that the susceptibility anisotropy of a given 
crystallite has a single-ion contribution, Axion(T) — 
Xion ~ Xtom an d a van Vleck contribution due to the 
matrix, Axb — X% ~ Xb '• ^ ne unspecified magnetic im- 
purities which contaminate our samples, and which are 
described by Eq. (0), are assumed to be isotropic, be- 
cause the the torque shown in Fig. 8 has no Curie-like 
upturn. Thus, we write the susceptibility anisotropy of 
a given sample as 



AxT P {x,T) = lx [xA Xion (T)+A XB (x) 



(9) 



where the temperature dependence of Axb is neglected. 
In the dilute limit, we approximate the van Vleck 
anisotropy by a linear expression Axb{x) — Axq + 



dAxBJx) 



dx 



to derive an 



x=0 



and use Eqs. (|8j) and 

approximate relation between the experimentally deter- 
mined quantities 






k l x (x-xi) 



-7"A X o, 



(10) 



where 7" is the alignement factor of the referent sam- 
ple. The experimental values obtained from Fig. 9 in- 



dicate that the ratio in Eq. (|Î0]) is indeed very nearly 
sample- independent, and given by 7^x0 — —2.4 • 
10~ 5 [emu /mol]. Using this value, and Eqs. (||) and (g), 
we can express the product 7^ A ^ i( _/T) in terms of ex- 
perimentally determined quantitiesJ23 



H\ X ion{T) = 



A x r P (x,T)-ll 71^x0 



X\ /C™ 



•i'l 



(11) 



Considering jfAxion(T) as a function of temperature 
one obtains the same function for all the samples. The 
value of 7" is obtained by fitting the universal curve (IO) 
by the theoretical result for the Coqblin-Schrieffer model, 
with the same parameters as in the previous section. This 
gives 7" = 0.55. 
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FIG. 10. The single-ion susceptibility anisotropy Ax*on(T) 
defined by Eq. ( |ll| ) plotted versus temperature for all the 
Ce I ,Lai_ I ,Cu2.05Si2 samples with x < 0.06. The full line is 
the scaling result, and the dashed line is the CF result. 

The single-ion anisotropy Axion(T), defined in the 
above described way, is plotted as a function of temper- 
ature in Fig. 10, together with the scaling result (full 
line) and the CF anisotropy (dashed line). The data 
show clearly that the magnetic anisotropy of 4/ 1 ions in 
the tetragonal Ce a; Lai_ :E Cu2.o5Si2 crystals can not be ex- 
plained by the CF theory, which considers only the ther- 
mal fluctuations among various CF states. The correct 
description is obtained by including the quantum fluctu- 
ations, and is provided here by the scaling solution of the 
Coqblin-Schrieffer model. 



C. Theoretical analysis 

The properties of Ce ions in dilute Ce :c Lai_ :c Cu2.o5Si2 
samples are modeled by the Coqblin-Schrieffer Hamilto- 
nian with CF splitting, 



H-k = Hqf + H o + Hj, 



(12) 



where Hcf describes the / states in the tetragonal CF 
(see the Appendix), Hq describes the conduction bând 
of width 2Dq, with a constant density of states p, and 
Hj defines the exchange-caupling between the 4/ 1 states 
and the bând electrons.EllE3 

The high-temperatuţe-properties are calculated by the 
poor man's scalingJljO which provides the renormal- 
ized coupling by reducing the conduction electron cutoff 
from Dq to -D, and simultaneously rescaling the coupling 
constant from Jo to J(D), so as to keep the low-energy 
excitations of the total system upchanged. Assuming the 
isotropic exchange, this leads toEâO 



pJexp(-- 



D 



k B T K 



D 



M 



(13) 



where J is the absolute value of the antiferromagnetic 
coupling constant, m — 2 and M — 4 are the degeneracies 
of the lower and the upper CF state, respectively, A is 
the CF splitting, and Tk is the Kondo temperature. The 
renormalized coupling constant J(D) defines an effective 
model with cutoff D. The Kondo temperature is given 
by Eq. (|1J) with J and D replaced by-the iniţial values 
Jo and D , respectively (for details seeE30). The scaling 
law (|Î3|) applies provided the renormalized coupling is 
small, i.e., pJ<S^l and D 3> Tk- 

The scaling trajectory described by Eq.(|i~3|) has two 
asymptotic regimes. The high-temperature asymptote, 
Jh{T), is obtained by neglecting the CF splitting, and is 
given by, 



p Jh exp 



pJh 



D 



N 



(14) 



where, TU is the high-temperature Kondo scale, defined 
by Eq.(|l4|) with J and D replacing Jh and D, respec- 
tively. 

The low-temperature asymptote Jl{D) is obtained by 
neglecting in Eq. ( |l3| ) the effects of the excited CF levels, 
i.e., by neglecting Tk and D with respect to A. This 
identifies Jl(D) as the scaling trajectory of an effective 
doublet with Kondo scale Tk, and the effective coupling, 



pJ L exp [ — 

pJh 



D 



(15) 



Eqs. ([[3]- |Î5|), show that the / level behaves at high tem- 
peratures as an effective sextet with Kondo scale Tj^ , and 
at low temperatures as an effective doublet with Kondo 
resonance Tk- The two Kondo scales satisfy the relation, 



(W 



H\m+M 
K ) 



{keT, 



K 



,A/ 



(16) 



The scaling result given by Eq.dÎ3) shows that we can 
remove the CF splitting from the asymptotic considera- 
tions at high and low temperatures, provided we renor- 
malize the effective coupling constant in such a way that 



the low-energy properties of the effective models and the 
Coqblin-Schrieffer model are the same. Note, the Tk 
of the effective doublet is much higher than the Kondo 
scale of the Coqblin-Schrieffer model calculated in the 
limit A 3> Do, with the iniţial condition Jl(Dq) = Jq. 

The anisotropic susceptibility tensor of the Coqblin- 
Scricffcr model with CF splitting is obtained by applying 
the linear response theory, and calculating the response 
functions by lowest-order pexiuxbation theory with renor- 
malized coupling constantsCJCJ. The effective Coqblin- 
Schrieffer model at a given temperaturc T is obtained by 
rcducing the cutoff from Do to D = Ak B T, and renor- 
malizing J to J{T) = J(T K /AT, A/Ak B T) according to 
Eq. (|l3|). Here, A is a numerical constant of the order of 
unity. Thus, we obtain the susceptibility 



xJ(T)=Xcf(T)[1-2 P J(T)}, 



(17) 



where Xcf(T) is the a component of the anisotropic 
CF susceptibility calculated in the Appcndix. The re- 
sult ( |l7| ) shows that the exchange interaction reduces 
the CF susceptibility, and the poor man's scaling gives 
the reduction factor as F(T,T K ,A) = [1 - 2pJ(T)]. 
Fitting the experimental results above 30 K with the 
renormalized susceptibility given by Eq.(|17|), where J(T) 
is given by Eq.(|13|), we obtain an implicit relation be- 
tween Tk and the cutoff constant A. The analysis of the 
Ce^Lai_ x Cu2.o5Si2 data above 30 K shows that we van 
write Tk — aA + (3, where a = 3 K and Ş — — 1 K. 
The physically acceptable range of the cutoffs is from 
k B T to 4k B T, which gives T K between 2 K and 11 K, 
and Tk between 60 and 110 K. However, the analysis of 
the low-temperature data (see below) leads to the values 
Tk = 8.5 and A = 3. The theoretical results obtained in 
such a way are shown in Figs. ||, [| |Î0[ [Tî], and [Î2| as a 
full line. 

The experimental data are described by the scaling 
theory down to 30 K. At lower temperatures there is a 
discrepancy, which is not really surprising, since J(T) in 
Eq.(§7|) grows rapidly for T < Tg and at about 30 K 
the coupling given by Eq.([l^) becomes too large for the 
renormalized perturbation expansion to be valid. Note, 
in the absence of the CF splitting the perturbation theory 
for the A^-fold degenerate / level breaks down at around 
T ~ Tjţ. The CF splitting reduces the degeneracy of the 
ground state and extends the validity of the perturbation 
expansion down to about T ~ ATk, where Tk is the 
Kondo temperature of the effective doublet, and Tk <C 

rpH 

1 k ■ 

Below 20 K, however, the coupling constants defined 

by Eqs.© and (|Î5|) satisfy J(T) ~ J L (T), such that we 

can consider the / state as an effective doublet. The av- 

erage moment of the CF ground state, calculated by the 

CF theory, is fx = 1.62/is, which is not too different from 

the free spin-1/2 moment, %/3/^s- Thus, to discuss the 

low-tempcrature data, for T < Tk we neglect the small 

anisotropy of the CF ground state and approximate the 

CF split multiplet by an effective doublet, and replace the 



Coqblin-Schrieffer model by an effective spin-1/2 Kondo 
model. We reiate the two models by setting the Kondo 
temperature of the spin-1/2 Kondo problem to Tk- In 
such way the parameter space of the low-temeperature 
model is complctly determined. The susceptibility of the 
Kondo model in the local moment (LM) regime, i.e., for 
Tk/2 < T < \5TKx-if given by the numerical renormal- 
ization group resultlia 



Xw = 



0.68C 
T + VÎTj 



(18) 



K 



where C = N A (gp B ) 2 S(S + l)/3fc s , i.e., C Q = 0.375 
emu/mol. Using \w an d Tk = 8.5 K we obtain the curve 
shown in Figs. 5 and 12. In the low-temperature LM 
regime, which sets in for diluite Ce a; Lai_ a ;Cu2.o5Si2 alloys 
between 4 K and 30 K, we find that the calculated sus- 
ceptibility is very close to the experimental values. Thus, 
by demanding that the low- and the high-temperature 
models have the same Kondo temperature we restricted 
the high-temperature cutoff constant to A = 3. 



D. Discussion of the experimental results and 
conclusion 

The 
anisotropic susceptibility of a single Ce ion is found by 
systematic measurements of dilute Ce x Lai^ K Cu2.o5Si2 
alloys on the Faraday balance and the torque magne- 
tometer, and by carefully subtracting the background. 
The experimental results are explained by the Coqblin- 
Schrieffer model with the CF doublet-quartet splitting 
A = 330 K and with the exchange scattering such that 
Tk = 8.5 K. The average value of the calculated suscep- 
tibility tensor, X j(T) = [2xf(T) + x}(î 1 )]/3, is shown 
in Figs. 5 and 6, together with the experimental data. 
The corresponding result for the susceptibility anisotropy 
Ax,j(T) is compared with the torque measurements in 
Fig. 



10 



Combining the Faraday balance and the torque 
data we find the single-impurity response along the prin- 
cipal axes, which is shown in Figjl^ for the x = 0.06 
sample, together with the scaling result (full lines) and 
the CF theory (dashed lines). The experimental data for 
other samples are about the same. 
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FIG. 11. The single-ion susceptibility for finite Ce concen- 
tration shown on a Curie- Weiss plot as a function of temper- 
ature for Ceo.o6Lao.94Cu2.o5Si2- Dashed line is trie CF result 
1/Xcf(T), and full line is the scaling result 1/xj( t )- The 
calculations are performed using Tk = 8.5 K, T^ — 100 K, 
and A = 330 K. 

The principal-axes susceptibilities shown in Fig. O fol- 
low between 100 K and 350 K an anisotropic Curie- Weiss 
law, 



C° 



Xion 



T -<d c 



(19) 



such that the slope of l/x? on ( T ) and 1 /Xcf( t ) is about 
the same. Thus, the high-temperature data can be dis- 
cussed in terms of an anisotropic local moment which is 
close to the CF value. However, the response of the / 
state is reduced with respect to the CF value due to the 
temperature dependent exchange coupling to the conduc- 
tion bând. The lowest order perturbation theory gives 
the reduction factor as F(T/T%) = 1 - 2pJ H (T/T§), 
where we assumed an isotropic exchange coupling. The 
Ce impurity behaves in this high-temperature LM regime 
as a sextet split by the tctragonal CF, and with the rel- 
evant Kondo scale T§ = 100 K. 

Below T ~ T^ , we observe the crossover from the 
high-temperature LM regime to an effective twofold de- 
generate low-temperature LM regime. Surprisingly, the 
behavior of the system at the crossover is still rather well 
described by scaling, and the single-ion response can be 
discussed in terms of the reduced CF susceptibility. 

At low temperatures, T < A Tk, the renormalized cou- 
pling becomes too large for the lowest-order renormalized 
perturbation expansion to be valid. Since we are not 
aware of any accurate theoretical results for the response 
of Coqblin-Schrieffer model with CF splitting that we can 
use close to Tk, we estimate the reduction factor F(T) di- 
rectly from the experiments, assuming that the exchange 
is isotropic and that the susceptibility retains the factor- 



ized form ([17]) down to lowest temperatures. The Fara- 
day balance data give F(T) — Xion(T)/xcF(T) and the 
torque data give F(T) — Axion(T)/AxcF(T), which we 
plot in Fig.[12], together with the high-temperature scal- 
ing result xj/xcf — [1 — 2pJ(T)] and the Wilson's result 
F w = Txw/Co- 
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FIG. 12. Temperature dependence of the reduction fac- 
tor of the susceptibility, F(T) = Xion(T)/xcF(T) for 
Ce :E Lai_a;Cu2.05Si2 with x < 0.09. The symbols are obtained 
from the susceptibility data, the dashed line is obtained from 
the anisotropy data and the solid line is the scaling result. The 
dashed-dotted line shows the reduction factor for the spin-1/2 
Kondo model. 

The reduction factors obtained from the average sus- 
ceptibility and the anisotropy data are very similar, and 
the high-temperature data are rather well described by 
the poor man's scaling. At low temperatures, where the 
crystal field calculations lead to an effective doublet with 
Xcf — N A pl s /3k B T and fi c g ^ 1.621ti S , the experi- 
mental reduction factor comes very close to the universal 
curve obtained for the isotropic spin-1/2 Kondo modcl.Eij 

We mention, for completeness, that the magnetic 
anomalies discussed here are aceompanied by the usual 
Kondo anomalies in the electric resistivity and the ther- 
moelectric power data.l3 AII the samples used for the sus- 
ceptibility measurements have a clear Kondo minimum in 
the bare resistivity, as well as a very large thermoelectric 
powe&.with two peaks, which is typical of CF split Kondo 
ions.E3 The low-temperature peak is at about 10 K, and 
the high temperature one at about 120 K, as expected in 
a system described by the Coqblin-Schrieffer model with 
the Kondo scales, T K ~ 8.5 K and Tg ~ 100 K. 

In summary, the magnetic susceptibility of a dilutc 
Ce ions embedded in the tetragonal mctallic host has 
been obtained by careful data analysis. The samples 
used in our studies have a negligibly small Ce-Ce inter- 
action. The changes in the matrix induced by the doping 



are also found to be very small. The behavior of the 
Ce s Lai_a;Cu2.o5Si2 compounds with less than 9 at% of 
Ce is well described by the Coqblin-Schrieffer model of a 
/ state comprising a ground state doublet and a pseudo- 
quartet split by A ~ 330 K. The quantum fiuctuations 
due to the exchange coupling between the / state and 
the conduction bând can be described by the poor man's 
scaling, which explains the high-temperature data and 
the crossover from the high-temperature LM regime to 
a twofold degenerate low-temperature LM regime, that 
takes place at about Tjţ ~ 100 K. The lowest-order per- 
turbation expansion, based on the scaling solution of the 
CF split model, breaks down below 30 K. However, at 
such low temperatures, the effect of the excited CF states 
is rather small, and the experimental results below 20 K 
can be described by the exact solution of the spin-1/2 
Kondo model with Kondo scale Tk — 8.5 K. This low- 
temperature Kondo scale determines completly the cutoff 
constant used in the high-temperature scaling. 
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F. Appendix 

A Ce 3+ ion in_a.tetragonal crystal field is described by 
the HamiltonianEa 



H C f — B 2 2 



BlOl 



B a Oa, 



(20) 



where B™ are the CF parameters and O™ are the Stevens 
operators, which are related to the angular momentum 
operator J and its components J x = \{J+ + </_), J y = 
y(J+ ~ J-)> and Jz> as follows: 



2 ° = 3J 2 -J(J + /), 



(21) 



Ol = 35Jf-30J(J+I) J2 + 25J2-6J(J+I)+3J 2 (J+/) 2 , 



and the energy eigenvalues e„ are 



ei = 10B° 2 + 60Bl + 12^5 



V 



T R 4 



e 2 = -2B° 2 - 180B 4 ° - U ^ 1 ^ Bl 



n 



(23) 



"1 IoulJ 4 

e 3 = -8B° + 120BI, 
where the mixing parameter r\ is defined by the equation 
t\ V 1 -»/ 2 _ B% + 20BI 



v/I - ^ 



V 



V5Bf 



(24) 



The neutron scattering datali.3 indicate a T7 ground 
state and a CF scheme shown in Fig. 13. 
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FIG. 13. The CF level for the doublet-quartet scheme used 
in the calculations. 

Sincc £2 ~ £3, we use here an approximate doublet- 
quartet scheme. 

The a component of the CF susceptibility of one mole. 
of isolated Ce 3+ ions is given by the van Vleck formulând 



a _ N A (gjfi B ) 2 ^ 
Xcf - v 2^ 



Ş J2 \{m\J a 



n)fe-P e ™ 



+ £ \{m\J a \n)\- 



-0e„ 



-0e n 



c n c m 



(25) 



oî 



;(4 



Jt 



For small CF this Hamiltonian is a perturbation to the 
sixfold degenerate 4/ , j — 5/2 atomic wave functions 
|±|) , |±f) , |±f) ; and it is easily diagonalized. The 
eigenvectors \n) are given by 



r 



(i) 



,(2) 



\±1)= V 
|±2) = - 
|±3) 



± 5 - 
2 



Vi 
1 



±- 






T 2 



=F- 



(22) 



where n,m £ {±1,±2,±3}, Na is the Avogadro num- 
ber, ^b the Bohr magneton, gj the Lande gyromagnetic 
factor, Z the partition function, [3 = l/fesT, and e n ,e m 
are given by Eq. (|23|). The matrix elements (m | J a | n) 
of the angular momentum J in a direction are taken be- 
tween the CF eigenstates (p2|). 

The n summation is performed over all the energy lev- 
els, while the m summation is performed for degenerate 
{(-n = £m) and for non-degenerate (e n ^ e m ) levels sep- 
arately. If the energies are measured relative to ei, we 
obtain 



Xcf 



2{l + 2e-l iA 



■Ş. 



2(§-V 
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r 



-f3A 



(26) 
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+21^- ^ 
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+ 64?y 2 (l-T] 2 



X& = ^&LP { [IO, 2 (1 ^ , 2 ) (1 + e-<") 



2(l + 2e-^ A )' 



"2 6 



-0A, 



877" 



(l-^)+5(l-2^) a 



-/3A 



A 



(27) 



The anisotropy Ax = x c ~ X ab i s 



Axcf 



N A {. gjl iB? {W - 1) 

2(l + 2e-^ A ) 



[IO77 2 (lArf - 17) e _/3A ] 13 



+ (1 V - 9) + (13 - 14r; 2 



1 - e^ A 



(28), 



In both, the low and the high temperature limits the 
CF susceptibility in a direction can be approximated by 
the Langevin formula, Xcf = ^A^a/3k B T 7 where fi a is 
the effective magnetic moment. At high temperatures, 
T ^ > A, the effective moment is isotropic, fi a b — fi c — 
9jVJU+T)Vb = 2. 535^ B . 

At low temperatures, T«A, the magnetic moment is 

determined by the Iy ground state doublet moment and 
it is anisotropic 



M. 



= V3gj 



-Arj' 



Hb, 



fi 



at = v / Î5.g,/?7Vl - V 2 Hb- 



(29) 



The curves shown in the text are calculated with pa- 
rameters A = 330 K, 77 = 0.816, and gj — 8/7. The cor- 
responding values for the CF parameters are B% = —1.01 
meV, Bl = 0.011 meV, and Bj = -0.50 meV. 

From those parameters we find fi c = 1.727 fi b and 
flab = 1.566/j.b. and the effective magnetic moment of 
the spherically averaged T 7 CF ground state p e ff = 

^2/3(^)2 + 1/3(Pc) 2 - 1.621flB. 
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samplc 


slopc k' x 


intcrcept l' x 


em-it 
mol 


li . X 
k\ 'o.Ol-x 


1-7 il)i- 

mol 


7i 


0.01 h 


0.566 


-1.4110 - 7 


- 


0.31 


0.02 v 


2.042 


2.73-10-° 


-2.67-10- 5 


0.56 


0.02 h 


0.901 


1.34-10- 6 


-2.97-10- 5 


0.24 


0.03 v 


-1.06 


1.71-10-° 


-2.42-10-° 


0.20 


0.03 h 


-2.78 


-4.45-10-° 


-2.40-10-° 


0.52 


0.06 v 


1.825 


3.61-10-° 


-2.37-10-° 


0.17 


0.06 h 


-4.170 


-8.46-10-° 


-2.43-10-° 


0.38 



TABLE I. The parameters from Fig. 9. 



12 



